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VERITAS observations of the Segue 1 dwarf spheroidal galaxy 
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Abstract: In the cosmological paradigm, Cold Dark Matter (DM) dominates the mass content of the Universe and is 
present at every scale. Candidates for DM include many extensions of the standard model, with a Weakly Interacting 
Massive Particle (WIMP) in the mass range from 50 GeV to greater than 10 TeV. The self-annihilation of WIMPs in 
astrophysical regions of high DM density can produce secondary particles including Very High Energy (VHE) gamma 
rays with energies up to the DM particle mass. The VERITAS array of Cherenkov telescopes, designed for the detection of 
VHE gamma rays in the 100 GeV- 10 TeV energy range, is an appropriate instrument for the detection of DM. Among the 
possible astrophysical targets, dwarf spheroidal galaxies (dSphs) of the Local Group are promising targets to search for 
the annihilation signature of DM due to their proximity and large DM content. We report here on extensive observations 
conducted by VERITAS on the nearby Segue 1 satellite galaxy, which is currently considered as one of the best dSphs 
for DM studies. The results are discussed in the framework of WIMP models, with a special emphasis on leptophilic DM 
models invoked to explain the recent cosmic-ray lepton anomalies. 



Keywords: dark matter; dwarf spheroidal galaxy; gamma-rays 



1 Introduction 

The compelling evidences for the presence of Dark Matter 
(DM) in the different structures of the Universe 1 1 1 have 
motivated numerous efforts to search for DM by means of 
astrophysical observations. If DM is made of Weakly In- 
teracting Massive Particles (WIMPs) annihilating into stan- 
dard model particles, indirect searches for DM annihila- 
tions with very high energy (VHE) 7-rays provide one of 
the best way to constrain the nature of the DM particle: 7- 
rays are free of any propagation effects on short distances 
(< 1 Mpc) and DM particle annihilation is predicted to give 
a unique 7-ray spectrum. Such searches are often con- 
ducted using pointed observations toward nearby DM over- 
densities, because the annihilation rate is proportional to 
the squared DM density. Popular targets include the Galac- 
tic Center 12 , 3j|4]j5], satellite galaxies of the Milky- Way 
(MW) I6ji7n8ii9l[l0][n][l2j, globular clusters and 
clusters of galaxies ifTSl . 

The dwarf spheroidal galaxies (dSphs) of the Local Group 
best meet the criteria for a clear and unambiguous detection 
of DM. They are gravitationally bound objects and con- 
tain up to (9(10'^) more mass in DM than in visible matter. 
As opposed to the Galactic Center, they are environments 
with a favorably low astrophysical 7-ray background. Nei- 
ther astrophysical 7-ray sources (supernova remnants, pul- 
sar wind nebulae, etc) nor gas acting as target material for 
cosmic rays have been observed in these systems lfT6l . Fur- 



thermore, their relative proximity and high galactic latitude 
make them the best astrophysical targets for high signal-to- 
noise detection. 

Segue 1 is an ultra faint dSph discovered in 2006 as an 
overdensity of resolved stars in the Sloan Digital Sky Sur- 
vey |,17J . It is located at a distance of 23 ± 2 kpc from the 
Sun at (Ra,Dec) = (10'"07'"03.2M6d04'25"), well above 
the Galactic plane. Because of its proximity to the Sagittar- 
ius stream, the nature of Segue 1 overdensity has recently 
been disputed, some people arguing that it was a tidally 
disrupting star cluster originally associated with the Sagit- 
tarius dSph 1 1 8 1 . However, a recent kinematic study of a 
larger member star sample (66 stars compared to the pre- 
vious 24 stars sample) has firmly confirmed that Segue 1 
is an ultra faint MW satellite galaxy I.19J . According to 
the study of its star kinematics. Segue 1 is the most DM- 
dominated dSph and is often highlighted as the best dSph 
target for indirect DM searches [,21. ,22. ,23] . 

2 Observations 

VERITAS [24] is an array of four 12 meter imaging at- 
mospheric Cherenkov telescopes (lACTs) located at the 
Fred Lawrence Whipple Observatory (FLWO) in southern 
Arizona. Each telescope is composed of a large miiTor 
area which reflects Cherenkov light induced by particle air 
showers on a camera placed in the focal plane. Each of the 
VERITAS cameras consists of 499 photomultiplier tubes 
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Figure 1: Significance map obtained after 7-ray selection and background subtraction. The significance is calculated 
according to the Li & Ma method |27|. The black cross indicates the target position. 



covering a 3.5° field of view. The large collection area 
(~10^ m^), in conjunction with the stereoscopic imaging 
of air showers, allow VERITAS to detect VHE gamma-rays 
between energies of 100 GeV and 30 TeV with an energy 
and angular resolution of ^15-25% and ~ 0.1°, respec- 
tively. A source with a 1% Crab Nebula flux can be de- 
tected by VERITAS in approximately 25 hours. 
Observations of the Segue 1 dSph were performed between 
January 2010 and May 2011, at a mean zenith angle of 
^ 20°. Only data taken under good weather conditions and 
with no major hardware problems were selected. The total 
exposure, after dead time correction, is nearly 48 hours. 
This is the biggest exposure reported so far by any lACT 
array on a dSph. The observations have been conducted us- 
ing the "wobble" pointing mode, where the camera center 
is offset by 0.5° from the target position. This allows for 
simultaneous background estimation and source observa- 
tion, thus reducing the systematic uncertainties in the back- 
ground determination. Four wobble directions were used 
(North, South, East, West) and were alternated from run to 
run. 

3 Data analysis and results 

Data reduction follows the methods described in [25|. Af- 
ter calibration of the camera gains, images recorded by 
each of the VERITAS telescope are characterized by a 
second moment analysis giving the Hillas parameters. A 
stereoscopic analysis combining each telescope image pa- 
rameters is used to reconstruct the 7-ray arrival direction 
and shower core position. Cuts on the Hillas parameters 
and the angle 9 between the target position and the recon- 
structed arrival direction (9^ < 0.015 deg^, defining the 
source region) are applied to reduce the cosmic-ray back- 
ground. The image distance from the center of the cam- 
era is also required to be less than 2° to avoid truncation 



effects. Furthermore, the integrated charge recorded in at 
least two telescopes is required to be > 90 photoelectrons 
which effectively sets the analysis energy thresholcO] to be 
170 GeV. 

After the 7-ray selection, the residual background was es- 
timated using the ring background technique [26|. Two 
circular regions, of radius 0.2° centered on the target po- 
sition, and of radius 0.3° centered on the bright star 77- 
Leonis (with apparent magnitude in the visible band V = 
3.5, and located 0.68° from the position of Segue 1), were 
excluded for the background determination. The resulting 
significance map is shown on figure[T] No significant 7-ray 
excess is found at any positions within the field of view. 
The large depletion area, with negative significances, cor- 
responds to 7]-Leonis. Given the absence of signal, one can 
derive an upper limit (UL) on the number of 7-rays in the 
source region. Because of the energy reconstruction bias in 
the low energy regime near the threshold, we define a mini- 
mum energy Emin — 300 GeV to safely compute ULs. Us- 
ing the Rolke method |28|, the UL at the 95% confidence 
level (CL) in the number of excess events in the source re- 
gion is: 

p^95%CL(-g ^ ggy^ ^ -^Q2.5 , (1) 

which gives an integral flux UL above 300 GeV of, assum- 
ing a power law spectrum with index F — 2.6: 

^95%CL(^g > 300 GeV) < 8 X 10"^^ cm"^ s^^ (2) 
This corresponds to 0.5% of the Crab Nebula integral flux. 

4 Dark Matter annihilation bounds 

The absence of signal at the position of Segue 1 can be used 
to derive constraints on various DM models. The 7-ray flux 

1. The energy threshold quoted here is taken to be the energy 
at which the differential detection rate of 7-rays from the Crab 
Nebula peaks. 
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from the annihilations of DM particles, of mass moM, in a 
spherical DM halo is given by a particle physics term times 
an astrophysics term: 



'^^ cm'^ s ^ around 1 TeV, 
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"dE 
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dE 



X 3{An), (3) 



where the astrophysical factor J(AJ7) is the squared DM 
density integrated along the line of sight (los) and over 
the solid angle Ail. The solid angle is given here by 
the size of the signal search region defined previously in 
our analysis, i.e. 9^ < 0.015 deg^. The particle physics 
term contains all the information about the DM particle: its 
mass mDM, its total velocity-weighted annihilation cross- 
section (crv) and the differential 7-ray spectrum from all 
final states weighted by their corresponding branching ra- 
tios, dN^/dE. 

The estimate of the astrophysical factor requires a model- 
ing of the Segue 1 DM profile. An Einasto profile \29 \ is 
used: 

PDM(r)-Pse-2"[('^A=)^'"-il, (4) 

with the scale density, the scale radius and the index n re- 
spectively being ps = 1.1 X 10^ Mq kpc~'^, Tg = 0.15 kpc 
and n = 3.3 |[T2 30 |. The value of the tidal radius changes 
the astrophysical factor by less than 10%. We adopt a value 
of 500 pc, which is the median truncation radius of the 
Via Lactea II simulation subhalos presenting similar char- 
acteristics to the Segue 1 DM halo lfT9ll20l . Having these 
parameters in hand, the value of the astrophysical factor 
within the solid angle subtended by our source region is 
3{Ail) = 7.7 X 10^* GeV"2 cm"^ sr. The systematic un- 
certainties quoted in the latest modeling of the Segue 1 DM 
distribution are less than an order of magnitude 121]. 
Once the astrophysical factor has been estimated, the UL 
in the number of 7-ray candidates in the source region can 
be translated into an UL on the total annihilation cross- 
section: 
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(5) 

where Tobs is the total observation time and Acb(E) is the 
effective area as a function of energy, zenith angle and the 
offset of the source from the pointing position. 

4.1 Classical constraints 

Independently of the DM model, we consider in this sec- 
tion "classical" annihilation channels, each time assum- 
ing a 100% branching ratio. Figure |2] shows our con- 
straints on ((Tv) as a function of the DM particle mass for 
three annihilation channels XX ~^ t~^t~, xX ~^ W''"W^ 
and XX ~^ bb. The 7-ray spectra have been simulated with 
the particle physics event generator PYTHIA v8.1 lOTl . 
The bb and r+r^ channels encompass all possible DM 
7-ray annihilation spectra in the VERITAS energy range, 
and give an idea of the uncertainties related to the particle 
physics model. For the W+W" channel, the 95% CL UL 



on ((Tv) is at the level of 8 x 10 
among the best reported so far with dSph observations. 
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Figure 2: Upper limits at the 95% CL on (fTv)o as a func- 
tion of the DM particle mass for different "classical" an- 
nihilation channels. The grey band area represents a range 
of generic values for the annihilation cross-section in the 
case of thermally produced DM. The red solid curve cor- 
responds to the Sommerfeld enhancement effect that could 
possibly arises for the W+W^ channel (see text for de- 
tails). 

We also consider the case where the annihilation cross- 
section can be boosted with the Sommerfeld enhancement 
effect. The Sommerfeld effect arises when the two DM par- 
ticles interact through an attractive potential, mediated by 
the exchange of a massive boson 1,32 J . The annihilation 
cross-section enhancement is particularly effective when 
the DM particle relative velocity is very small. Depend- 
ing on the mass and the coupling a of the exchanged 
boson, the Sommerfeld enhancement can exhibit a serie of 
resonances for specific values of the DM particle mass, giv- 
ing very large boost factors up to 10^. In such scenarios, 
the annihilation cross-section is given by: 



(crv) = S X (crv)o, 



(6) 



where S is the Sommerfeld boost factor averaged over the 
DM particle velocity distributior0 and depends on the DM 
particle mass, (ctv)o is the thermal WIMP annihilation 
cross-section before freeze-out. The Sommerfeld enhance- 
ment is of particular interest for cold DM halos like dSphs, 
where the mean star velocity dispersion can reach 10~^ c. 
Figure |2] shows the limits on (fTv)o, accounting for the 
corresponding Sommerfeld boost in the case of DM an- 
nihilating into a pair of W+ bosons when a Z° boson 
{mza ^ 90 GeV, a ^ 1/ 30) is exchanged. The exclusion 
curve exhibits resonances, which correspond to the first res- 
onances in the Sommerfeld enhancement for DM particle 
masses of 4.5 and 17 TeV, respectively. VERITAS severely 
excludes these first resonances. 



2. In the case of Segue 1, we assume that the relative velocity 
distribution is Maxwellian (i.e. the DM gaz is thermalized and in 
equilibrium), with a width of the order of the mean star velocity 
dispersion. See |21 1 for further details. 
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4.2 Models with a new force in the dark sector 

The last section is dedicated to a class of "leptophilic" 
DM models invoked to explain the recent cosmic -ray lep- 
ton anomalies measured by ATIC 133] and PAMELA ll34l . 
In any DM interpretation, the cosmic -ray lepton excesses 
require a DM particle mostly annihilating into leptons 
and with a high annihilation cross-section (compared to 
((Tv)o 3 X 10^^^ cm'^ s^^). Models which introduce a 
new force in the dark sector have recently been proposed 
to fit these requirements ll35l . The new force is carried by 
a light scalar field cj) predominantly decaying into leptons 
and with a mass and coupling to standard model particles 
chosen to prevent the overproduction of antiprotons. These 
models naturally comprise a Sommerfeld enhancement to 
the annihilation cross-section. As for the classic W+W^ 
channel, the boost also exhibits a serie of resonances, but 
for different DM particle masses because the mass of the 
exchanged boson and its coupling constant are different. 
Figure |3] shows the constraints on such models, for the two 
leptonic annihilation modes XX (f><P e+e^e+e~ and 
XX —^4"P~^ M^/i^M^/^^, and for a scalar particle (p of 
mass 250 MeV. VERITAS excludes most of the resonances, 
especially for the e^/e" channel. This result holds and 
strengthens for larger cj) masses, thus disfavoring this class 
of models. 
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Figure 3: 95% CL upper limits on {<jv)o as a function of 
the DM particle mass for lepton annihilation channels, in 
the framework of particle physics models with a new force 
in the dark sector The greay band area in the same as in 
figure |2] The dashed curves correspond to limits without 
any enhancement to (crv)o. The solid curves includes the 
Sommerfeld effect expected within these models. 
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